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Abstract 



The work functions (WF) of single-walled carbon nanotubes and bundles 
are studied using first principles methods. For individual metallic tubes, the 
WF is independent of the chirality and increase slightly with tube diame- 
ter. For semiconducting tubes, the WF (as defined by the HOMO energy) 
decreases rapidly. The WF of nanotube bundles (~ 5 eV) shows no clear de- 
pendence on the tube size and chirality, slightly higher than individual tubes. 
Calculations on finite tubes show no substantial difference in the tube end and 
the side wall. Upon alkali-metal intercalation, the WF decreases dramatically 
and the electronic states near the Fermi level are significantly modified. The 

metallic and semiconducting nanotubes bundles become indistinguishable. 
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The work function is one of the critical quantities in understanding the field emission 
properties of carbon nanotub esU. Although the work function can be estimated from the 
field-emission spectra based on Fowler-Nordheim mo delil, the results are not reliable due 
to the uncertainty of the local geometry of nanotubes@. Other experimental measurements 
on work functions of both single-walled nanotubes (SWNTs) and multi-walled nanotubes 
(MWNTs) included ultraviolet photoemission spectroscopy and transmission elec- 

tron microscopylii. It was found that the work functions of MWNTs are about 0.1-0.2 eV 
lower than that of the graphitei'00lli, while the SWNT bundles have slightly higher work 
functionJii. Upon Cs intercalation the work functions of carbon nanotubes are reduced 
dramaticallyi'0, which leads to a significant enhancement in field emissioni. Up to now, 
theoretical works were only limited to a few nanotubes with finite-leng thsM. In this letter, 
we report results of first principles calculations on the work functions of individual carbon 
nanotubes and nanotube bundles. The effect of tube diameter and chirality, and the differ- 
ence between the tube end and the side wall are investigated. The alkali-metal (K, Rb, Cs) 
intercalated carbon nanotube bundles are also discussed. 

The work function of a bulk metal is related to its Fermi energy Ep by WF = — Ep, 
where is the electrostatic potential caused by "spilling out" of electron density at the metal 
surface@^. For those metals with low electron density such as K, Rb, Cs, it is known that 
is much smaller than Ep §0. For carbon nanotubes, the conduction electron density is much 
smaller than that of K or cM. In this work, we approximate the WF by the Fermi energy 
To determine the Fermi level of carbon nanotubes with respect to the vacuum level, 
we perform all electron LCAO calculations based on the density functional DMol packaged. 
The density functional is treated by the generalized gradient approximation (GGA)i with 
the exchange-correlation potential parameterized by Wang and Perdew (PW91)0. For the 
infinite nanotubes, one- dimensional (1-D) periodic boundary condition is applied along the 
tube axis. The Brillouin zone is sampled by large sets of Monkhorst-Pack k meshesii (Along 
the tube axis, 40 k points are used for standard calculations and 160 points for accurate 
electronic density of states). Benchmark calculations are carried out on several alkali metals 



solids and the graphene sheet. The ionization potential of a Ceo cluster is also calculated 
from the total energy difference between the neutral and the charged cluster. As shown in 
Table I, the overall agreement between theory and experiment si'iSil for different systems 
is reasonable. 
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TABLES 

Table I. Work function of bulk alkali metals (K, Rb, Cs)il, graphene sheetlH, nanotube rope 
(see Table II and discussions in text)iJ^, and the ionization potentials of Geo cluster@. GGA 
denotes present GGA calculations, Exp. are the experimental WF valuesil^lBil. All units 
are in eV. 





K 


Rb 


Cs 


graphene 


nanorope 


Ceo 


WF (GGA) 


2.69 


2.42 


2.31 


4.8 


~ 5.0 


7.87 


WF (Exp.) 


2.30 


2.16 


2.14 


4.8i, 4.6i 


5.05!, 4.8i 


7.61 



We first address the work functions of the individual metallic SWNTs of infinite length, 
with diameter ranging from 5.7A to 16.3 A . Both the armchair (n,n) and the zigzag (n,0) 
{n = 3m) chirality are considered. The calculated WFs range from 4.63 to 4.77 eV, given in 
Fig.l and Table II. Gareful examinations show that the WF of the metallic tubes increases 
slightly with 1/D (Fig.l). Extrapolation towards the D— > 00 limit gives a WFqo of 4.83 eV, 
close to the calculated value for the graphene sheet (4.8 eV). It is worthy to note that the 
WFs of both the armchair and the zigzag metallic tubes fit the same linear dependence, 
indicating that the WFs of metallic nanotubes are independent of chirality. 
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FIG. 1. Work functions (eV) of individual metallic and semiconducting SWNTs vs. the in- 
verse tube diameter 1/D (squares: armchair SWNTs; dots: zigzag SWNTs). Linear extrapolation 
towards the D— > oo limit yields WFqo of 4.84 eV for metallic tubes and 4.73 for semiconducting 
ones, close to the calculated value (4.8 eV) for graphene sheet. 

For the purpose of comparison, we define the WF of semiconducting tube as the highest 
occupied molecular orbital (HOMO) energy. The WFs of the semiconducting tube are sub- 
stantially higher than the metallic ones (Fig.l). The calculated WF decreases linearly with 
1/D and approaches an extrapolation limit 4.73 eV at D— > oo. The strong diameter depen- 
dence come from the well-known gap decrease with tube diameter in semiconducting tubes0. 
There is no direct experiment measurement on the work functions of individual SWNTs. 
But recent experiments on the MWNT tips suggest that the WFs of semiconducting tubes 
(~ 5.6 eV) are high than those of metallic ones (~ 4.6 - 4.8 eV)0. 

To explore the difference in work function between the side wall and the tube end, we 
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carry out model calculations on several finite (5,5) tubes with length up to 21.8 A. As shown 
in Fig.2, one end of the nanotube is capped by a half-Ceo while the other end is H-terminated. 
The calculated HOMO energy decreases with the tube length I and approaches the infinite 
limit rapidly (4.87 eV, 4.63 eV, 4.55 eV for I = 3.8A, 11.7A, 21.8A respectively). In Fig.3, we 
plot the HOMO wave function distributions for the longest tube studied (160 carbon atoms 
on side wall plus 30 carbon atom on the end). The HOMO are clearly delocalized over the 
whole tube. The wave function at the tube end is comparable to that on the side wall. 
Similar distributions can be found for the molecular orbitals just below the HOMO level. 
Therefore, we expect that the work function at the tube end should not be much different 
from that of the sidewall. For the finite size nanotubes investigated, we find interesting 
oscillations in both the HOMO energy and the wave function distribution with the tube 
length. The oscillation period is three times that of the unit cells length along the tube 
axis, corresponding to the Fermi wavelength of armchair tubesS. This oscillation behavior 
implies that care should be taken when extrapolating the finite size results to the infinite 
limitM. Our detailed investigations on the electronic properties of finite nanotubes with 
different lengths and ends will be presented elsewhere. 




FIG. 2. Distribution of HOMO wave function (slice through the middle of nanotube) in a finite 
length (21.8 A ) tube with one closed end. The red (blue) color denotes positive (negative) value of 
the wave function. The delocalization of HOMO is clearly shown. The oscillations of wave function 
along the tube axis with a period of 3 unit cells correspond to the Fermi wavelength of armchair 
tubes. 



Table II. Work function (eV) of individual metallic tubes and ropes (tube bundle) with 
various diameters and chirality (m,n) obtained from GGA calculations. 



(m,n) 


(5,5) 


(6,6) 


(12,0) 


(8,8) 


(15,0) 


(10,10) 


(12,12) 


tube 


4.68 


4.71 


4.73 


4.74 


4.73 


4.76 


4.77 


rope 


5.08 


5.07 


5.05 


5.00 


4.98 


5.01 


4.94 
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We now discuss the work functions of the nanotube bundles. The bundles are modeled by 
two dimensional hexagonal lattices of uniform SWNTsH. For all the metallic tube bundles 
studied, the calculated work functions are around 5 eV (see Table II), slightly higher than 
individual metallic nanotubes and the graphene. The increase of WFs in tube bundles can 
be understood by the tube-tube interactional. From the experiments based on UPS, the 
work functions of SWNT bundles are determined to be 4.8 evi or 5.05 evU. Our theoretical 
results agree well with the available data on nanotube bundles. 

In addition to pristine materials, the electronic properties of carbon nanotubes can be 
efficiently controlled by alkali-metal intercalations. Recently, nanoelectronics devices have 
been constructed on the basis of alkali-metal doped carbon nanotube^. Thus, it is impor- 
tant to investigate the effect of alkali-metal intercalations. We have carried out systematical 
calculations on the effect of intercalations in both metallic (10,10) and semiconducting (17,0) 
tubes bundles. Intercalation density up to Kq.iC (Rb, Cs) (close to the saturation density 
in nanotube bundles^ and the graphite@) are studied. Both the inside of SWNTs and the 
interstitial sites are explored. The initial configurations of alkali-metal atoms are chosen to 
maximize the ion-ion distance^. Structural relaxations find no significant change from the 
initial configurations. For a given intercalation concentration, the WF is insensitive to the 
detailed configurations of the intercalated atoms. The lattice constants of two-dimensional 
hexagonal lattices are expanded. For example, we find that the intercalations of K atoms into 
the interstial sites of (10,10) tube bundle can induce about ~2A lattice expansion, which 
is comparable to expansion of 1.95 A in K intercalated graphite^ and 1.85 A in the HNO-, 
intercalated SWNT bundled^. In Fig. 3, we present the calculated work functions of interca- 
lated tube bundles as functions of the metal concentration for various metals. In general, the 
work function dramatically decreases after small amount of alkali-metal intercalations. The 
reduction becomes much slower at higher intercalation density. Furthermore, there is almost 
no difference between the WFs of intercalated (10,10) and (17,0) tube bundles. Thus, the 
metallic and semiconducting tube bundles become indistinguishable. The WF for Cs-doped 
nanotube bundles is slightly lower than that of the Rb- or K- doped systems. This trend from 
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K to Cs is consistent with the WFs for bulk metal that (WF{K) > WF{Rb) > WF{Cs), 
see Table I). 




FIG. 3. Work functions (eV) of alkali-metal intercalated carbon nanotube bundles vs. the 
intercalation density {x in M^C, M denotes metal). Initially, work function decreases dramatically 
with intercalation density. The reduction of WF becomes much slower at higher density. 

Intuitively, the reduction of WFs can be understood by the charge transfer from metal 
to carbon nanotube, which shifts the Ep of conduction band towards the vacuum level. 
Experimentally, the charge transfer from metal atoms to the nanotube have been confirmed 
by the resistance@'@ and Raman spectrumil. However, our present results show that in- 
teraction between K (or Rb, Cs) and nanotube cannot be simply described by a rigid band 
model with charge transfer. Shown in Fig.4 is the electronic density states for pristine and 
K-intercalated (10,10) SWNT bundles. Similar to the case of Li intercalationll, the valence 
bands are almost not affected by K intercalations. In contrast, the conduction bands are 
significantly modified by potassium-carbon interaction. New peaks associated with alkali- 
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metal atoms are found in the conduction bands. The density of states near the Fermi level 
is significant enhanced by the contributions of alkali metals. Our thoeretical results are 
supported by recent experiments on the optical properties of the K or Cs intercalated nan- 
otube bundles^. In addition, we have also calculated the density of states for intercalated 
semiconducting bundles. Similar to the case of work function, the metallic and semicon- 
ducting nanotubes bundles become indistinguishable. This is consistent with recent NMR 
measurements on K-intercalated SWNT bundlesil. 
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FIG. 4. Electronic density of states of pristine (dotted line) and K-intercalated (solid line) 
(10,10) SWNT bundles (K0.05C). The valence bands of nanotube are almost not affected by K 
intercalations, while the conduction bands are significantly modified by the potassium-carbon in- 
teractions. 

In summary, we have performed first principles calculations on the work functions of 
pristine and intercalated SWNT nanotube and bundles. The WFs of metallic nanotubes 
are weakly depend on the tube size and comparable to the graphene sheet. The work 



10 



functions of all the metallic tube bundles are around 5 eV, slightly higher than those of 
individual tubes. Finite size tube calculations reveal the delocahzed HOMO state, implying 
no substantial difference between the WF of the tube end and the side wall. Upon alkali- 
metal intercalation, the WFs of tube bundles decrease dramatically and the electronic states 
near Fermi level are significantly modified. The metallic and semiconducting nanotubes 
bundles become indistinguishable. This work is supported by the U.S. Army Research 
Office Grant No. DAAG55-98-1-0298, the Office of Naval Research Grant No. N00014-98- 
1-0597 and NASA Ames Research Center. The authors thank Prof. O.Zhou, Prof. Y.Wu 
for helpful discussions. We acknowledge computational support from the North Carolina 
Supercomputer Center. 

*E-mail address: zhaoj@physics.unc.edu 
^E-mail address: jpl@physics.unc.edu 
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